11 Perturbation Theory and Feynman Diagrams 



Wc now turn our attention to interacting quantum field theories. All of the 
results that we will derive in this section apply equally to both relativistic and 
non-relativistic theories with only minor changes. Here we will use the path 
integrals approach we developed in previous chapters. 

The properties of any field theory can be understood if the iV-point Green 
functions are known 



GNixi,...,XN) = {0\Tq^ixi)...c^ixN)\0} (1) 

Much of what we will do below can be adapted to any field theory of interest. 
We will discuss in detail the simplest case, the relativistic self-interacting scalar 
field theory. It is straightforward to generalize this to other theories of interest. 
We will only give a summary of results for the other cases. 



11.1 The Generating Functional in Perturbation Theory 

The A^-point function of a scalar field theory, 

Gn{xi,...,xn) = {0\T(j){xi)...q){xN)\0), (2) 
can be computed from the generating functional Z[J] 



Z[J] = {0\Te J |0) (3) 

In D = d + 1-dimensional Minkowski space-time Z[J] is given by the path 
integral 

r iS[(j)] + i I d^xJ{x)(l){x) 
Z[J] = V(f>e J (4) 



where the action S[(j)] is the action for a relativistic scalar field. The A^-point 
function, Eq.(l), is obtained by functional differentiation, i.e., 

1 (5^ 

Gn{xi,...,xn) = , , Yj-, ^Z[J] (5) 

Z[J\ dJ(xi) . . .dJ[XN) ./=o 

Similarly, the Feynman propagator Gf{xi —X2), which is essentially the 2-point 
function, is given by 

Gpixi - X2) = -^(O|T0(xl)0(a.2)|O) = I f Z[J] (6) 

Z[J\ Oj{xi)dJ{X2) ,7=0 

Thus, all we need to find is to compute Z[J]. 

We will derive an expression for Z[J] in the simplest theory, the relativistic 
real scalar field with a 0^ interaction, but the methods are very general. We 
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will work in Euclidean space-time {i.e., in imaginary time) where the generating 
function takes the form 



Z[J] 



V(j) e 



-S[(j)] + j d^x 



J{x)(j){x) 



where S[(j)\ now is 

S[(j>] = j d'^x 
In the Euclidean theory the A''-point functions are 

Gn{xi,. . . ,Xn) = {<l>{xi) . . .(j){xN)) 



Z[J] 6J{xi)...6J{xn) 



Z[J] 



j=o 



(7) 



(8) 



(9) 



Let us denote by Zo[J] the generating action for the free scalar field, with action 
5o [(/)]. Then 



Zo[J] 



I 



V(f) e 



d^xJ{x)(t>{x) 



= [Det (-a^ + 



/2 1^ d^x d''yJ{x)Goix - y)J{y) 



(10) 



where is the Laplacian operator in £)-dimensional Euclidean space, and 
Go{x — y) is the free field Euclidean propagator [i.e., the Green function) 



Go{x -y) = {(t>{x)(t){y))o = {x\- 



1 



-d"^ + m'- 



■\y) 



(11) 



where the sub-index label denotes a free field expectation value. 

We can write the full generating function Z[J] in terms of the free field 
generating function Zo [J] by noting that the interaction part of the action con- 
tributes with a weight of the path-integral that, upon expanding in powers of 
the coupling constant A takes the form 



-j d''x^^\x) 



= E 



(-ir 



n=0 



d^xi... / d^a;„ <p\xi)...<t>\x^) 



(12) 
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Hence, upon an expansion in powers of A, the generating function Z[J] is 



Z[J] = V<f)e 



= E 



(-1)" /A\" 



„ n! V 4! 



J d^x, ...j d^x„ j Vcp ct>\x,) . . . ct>\x„) e J 



= E 



(-1)" f XY f n f n S'^ 



= e 



4! y " sj{xr z,[j] 



(13) 



where the operator of the last Une is defined by its power series expansion. We 
see that this amounts to the formal replacement 

Sint[<f>] ^ Sint (jj^ (14) 

This expression allows us to write the generating function of the full theory Z[J] 
in terms of Zo [J] , the generating function of the free field theory, 

Z[J] = e (^■^) Zo[J] (15) 

Notice that this expression holds for any theory, not just a (j)'^ interaction. This 
result is the starting point of perturbation theory. 

Before wc embark on explicit calculations in perturbation theory it is worth- 
while to see what assumptions we have made along the way. We assumed, (a) 
that the fields obey Bose commutation relations, and (b) that the vacuum (or 
ground state) is non-degenerate. 

The restriction to Bose statistics was made at the level of the path integral 
of the scalar field. This approach is, however, of general validity, and it also 
applies to theories with Fermi fields, path integrals over Grassmann fields. As 
we saw before, in all cases the generating functional yields vacuum (or ground 
state) expectation values of time ordered products of fields. 

The restriction to a non-degenerate vacuum state has more subtle physical 
consequences. We have mentioned that, in a number of cases, the vacuum may 
be degenerate if a global symmetry is spontaneously broken. Here, the thermo- 
dynamic limit plays an essential role. For example, if the vacuum is doubly 
degenerate, we can do perturbation theory on one of the two vacuum states. 
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If they arc related by a global symmetry, the mimber of orders in perturba- 
tion theory which are necessary to have a mixing with its degenerate partner is 
approximately equal to the total number of degrees of freedom. Thus, in the 
thermodynamic limit, they will not mix unless the vacuum is unstable. Such 
an instability usually shows up in the form of infrared divergent contributions 
in the perturbation expansion. This is not a sickness of the expansion but the 
consequence of having an inadequate starting point for the ground state! In 
general, the "safe procedure" to deal with degeneracies that are due to symme- 
try is to add an additional term (like a source) to the Lagrangian that breaks 
the symmetry and to do all the calculations in the presence of such symmetry 
breaking term. This term should be removed only after the thermodynamic 
limit is taken. 

The case of gauge symmetries has other and important subtleties. In the 
case of Maxwell's Electrodynamics we saw that the ground state is locally gauge 
invariant and that, as a result, it is unique. It turns out that this is a generic 
feature of theories that are locally gauge invariant for all symmetry groups and in 
all dimensions. There is a very powerful theorem (due to S. Elitzur, and which 
we will discuss later) which states that not only the ground state of theories 
with gauge invariancc is unique and invariant, but that this restriction extends 
to the entire spectrum of the system. Thus only gauge invariant operators have 
non-zero expectation values and only such operators can generate the physical 
states. To put it differently, a local symmetry cannot be broken spontaneously 
even in the thermodynamic limit. The physical reason behind this statement is 
that if the symmetry is local it takes only a finite order in perturbation theory 
to mix all symmetry related states. 

Thus whatever may happen at the boundaries of the system has no con- 
sequence on what happens in the interior and the thermodynamic limit does 
not play a role any larger. This is why we can fix the gauge and remove the 
enormous redundancy of the description of the states. Nevertheless we have 
to be very careful about two issues. Firstly, the gauge fixing procedure must 
select one and only one state from each gauge class. Secondly, the perturba- 
tion theory is based on the propagator of the gauge fields i(0|TA^(x)Ay (a;')|0) 
which is not gauge invariant and, unless a gauge is fixed, it is zero. If a gauge 
is fixed, this propagator has contributions that depend on the choice of gauge. 
But the poles of this propagator do not depend on the choice of gauge since they 
describe physical excitations, e.g., photons. Furthermore although the propaga- 
tor is gauge dependent it will only appear in combination with matter currents 
which are conserved. Thus, the gauge-dependent terms of the propagator do 
not contribute to physical processes. 

Except for these caveats, we can now proceed to do perturbation theory for 
all field theories of interest. 

11.2 Perturbative Expansion for the Two-Point Function 

Let us discuss the perturbative computation of the two-point function in <j>'^ 
field theory in £)-dimensional Euclidean space-time. Recall that under a Wick 
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rotation, the analytic continuation to imaginary time ixo xd, the two-point 
function in D-dimensional Minkowski space-time, (0|T(/)(2;i)(/>(a;2)|0) maps onto 
the two-field correlation function of £>-dimensional Euclidean space-time, 



(o|r0(xi)0(x2)|o) 



{(t){xi)(j){x2)) 



(16) 



Let us formally write the two point function G^'^^xi — x^) as a power series 
in the coupling constant A, 



G^^\X^-X2) 



y:^_g^\x,-x,) 



n=0 



However, using the generating functional Z[J] we can write 

1 ^2 



G(2)(xi-a;2) = 



Z[J] SJ{xi)6J{x2) 



Z[J] 



7=0 



where 



(17) 



(18) 



(19) 



Z[J]=e V^-^/ Zo[J] 

Hence, the two-point function can be expressed as a ratio of two series ex- 
pansions in powers of the coupling constant. The numerator is given by the 

expansion of Z[J] 



j=o 



5J{xi)5J{x2) 



Z[J] 



j=o 



E 

n=0 



(-1)' 



5J{xi)5J{x2)6J{yiy' 5J{ynY 



Zo[J] 



j=o 
(20) 



and the denominator by the expansion of Z[0], which leads to a similar ex- 
pression but without a contribution due to the external legs, corresponding to 
the functional derivatives with respect to the source at the external points Xi 
and X2- The equivalent expression in Minkowski space- time is obtained by the 
replacement, 

-A < — > iX (21) 

at every order in the expansion. 

We will now look at the form of the first few terms of the expansion of the 
two-point function in perturbation theory. 



11.2.1 Zeroth Order in A. 

To zeroth order in A {i.e., 0(A°)), the numerator reduces to 



SJ{xi) 5J{X2) ^'"^^ J=o 6J{xi) SJ{x2) ^"'"^^ 

= Go{xi - X2) + 0{X) 



.7=0 



0(A) 



(22) 
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while the denominator is simply equal to one 

Z[0] = Zo[0] + 0{X) = 1 + 0(A) 

Hence, 

G'-^\xi - X2) = Go{xi - X2) + 0(A) 
11.2.2 First Order in A. 

To first order in A, the denominator Z[0] is given by 



(23) 
(24) 



Z[0] = 1 



-(- 

1! W! 



dy- 



5^ 



-Zo[J]|j=o + 0(A2 



(25) 



The expression inside the integrand can be calculated from the Taylor expansion 
of Z[)[J] in powers of J{x). To find a non-zero contribution wc need to bring 
down from the exponent enough factors of J so that they can be cancelled by 
the functional derivatives. Since the argument of the exponential factor in Zo[>^] 
is bilinear in J{x), 

1 



Zo[J] = [Det(-a2 + m2 



-1/2 ^ 2 



d^x / d^yJix)Goix-y)Jiy) 



(26) 



only an even number of derivatives in J{x) can be cancelled to a given order. 
In particular, to first order in A, we have to cancel four derivatives. This means 
that we need to expand the exponential in Zq [J] to second order in its argument 
to obtain the only non- vanishing contribution to first order in A to Z[J] at J = 0, 

Z[Q] = 1 + 



til 
1! 



A 
4! 

+0(A2 



5J{xY 2! 



d^'yi / d''y2J{yi)Go{yi - y2)J{y2. 



(27) 



The derivatives yield a set of 5-functions 



SJ{x) 



:[J{yi)---J{y4)] 



J=0 



p j=i 



(28) 



where P runs over the 4! permutations of the four arguments yi,y2, ya and y^. 
We can now write the first order correction to Z[0], in the form 

Z[0] = 1 + 

tt) (i) i (^) V "^"''^ / '^"^^ ' ■ ■ "^""y'^'^y' - y'^^'^y' ^(yp^ - ^' 



+0(A2) 



1 + 5 y d""^ Go{x, x) Go{x, x) + 0(A2) 



(29) 
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where 

4! 



^ = 2!2^ - 3 (30) 



'4 



Figure 1: A vertex and its four contractions. 

It is useful to introduce a picture or diagram to represent this contribution. 
Let us mark four points yi, ... ,1/4 and an additional point at x (which we will 
call a vertex) with four legs coming out of it. Let us join yi and 2/2 by a line 
and j/3 with by another line. To each line we assign a factor of Go{yi — 2/2) 
and Go(j/3 — 2/4) respectively, as in the figure: 

X y 

Goix - y) = 



Next, because of the ^-functions, we have to identify each of the points 
j/i , . . . , 2/4 with each one of the legs attached to y in all possible ways (as shown 
in figure 1) 

The result has to be integrated over all values of the the coordinates and of 
X. The result is 

Z[0] = 1 - (^) / (Go(x, x)f + 0(A2) (31) 

Physically, the first order contribution represents corrections to the ground state 
energy due to vacuum fluctuations. This expression can be represented more 
simply by the Feynman diagram shown in Fig. 2. 

Here, and below, we denote by a full line the bare propagator Go{x — y). 

Let us now compute the first order corrections to - — - — — — - — - Z\J] 

oJ{xi)oJ(x2) J=o 
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Figure 2: Vacuum fluctuations: first order correction to the denominator, Z[0]. 



We obtain 
(-1) A 



1! Alj"^ y 5J{xi)5.J{x2) 5J{yY ^^'^^ 



.1=0 



The non-vanishing contributions are obtained by matching the derivatives in 
Eq.(32) with an equal number of powers of J. We see that we have six factors 
of the source J at points zi, . . . ,zq and six derivatives, one at xi and at X2, 
and four at y. To match derivatives with powers amounts to find all possible 
pair-wise contractions of these two sets of points. Once the delta functions 
have acted we are left with just one integral over the position y of the internal 
vertex. Hence, the result amounts to finding all possible contractions among 
the external legs at Xi and X2, with each other and/or with the internal vertex 
at y. Notice that for each contraction we get a power of the bare (unperturbed) 
propagator Go- 

The only non- vanishing terms resulting from this process are represented by 
the Feynman diagrams of Fig. 3. 

OO 

X, y y + 

^ X2 . Xl X2 Xl 



7=0 

(32) 



X 



(a) 




Figure 3: First order contribution to the two-point function. 



The first contribution. Fig. 3(a), is the product of the bare propagator Go{xi — 
X2) between the external points and the first order correction of the vacuum di- 
agrams: 



(a) = Go{xi - X2) X 



J d'^y (Go(y,2/) 



(33) 
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The second term, the "tadpole" diagram of Fig.3(b), is given by the expression 



(b) = -[^) S I d^y Go(xi,2/)Go(2/,y)G)(2/,X2) (34) 
where the multipHcity factor S is 

5 = 4x3 (35) 

It counts the number of ways (12) of contracting the external points to the 
internal vertex: there are four different ways of contracting one external point to 
one of the four lines attached to the internal vertex at y, and three different ways 
of contracting the remaining external point to the one of the three remaining 
lines of the internal vertex. There is only one way to contract the two leftover 
internal lines attached to the vertex at y. 

By collecting terms we get the result shown in Fig. 4. 

To first order in A, the expansion of the two-point function can be written 
in the form shown in Fig. 5. 

At least at this order in perturbation theory, a number of diagrams which 
contribute to the expansion of the numerator get exactly cancelled by the expan- 
sion of the denominator, the vacuum diagrams. The diagrams that get cancelled 
are unlinked in the sense that one can split a diagram in two by drawing a line 
the does not cut any of the propagator lines. These diagrams contain a factor 
consisting of terms of the vacuum diagrams. We will see shortly that this is a 
feature of this expansion to all orders in A. 

Thus, to first order in A, the two-point function is given by (see figure) 

G^^\x^,X2) = G^^\x^,X2) -^Jd^'y G^^\x,,y)G^^\y,y)G^^\y,X2) + 0{X^) 

(36) 

as shown in Fig. 6. 

11.3 Cancellation of the vacuum diagrams 

The cancellation of the vacuum diagrams is a general feature of perturbation 
theory. Let us reexamine this issue in more general terms. We will give the argu- 
ments for the case of the two-point function, but they are trivial to generalize to 
any iV-point function. This feature also holds in all theories, relativistic or not, 

1 + oo 

Figure 4: Feynman diagrams for the two-point function to first order in A. 
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( 1+ OO +oi^')) ( + w +Q(A^)) 

G(2) = 



( 1+00+0(A2)) 



Figure 5: Factorization of Feynman diagrams for the two-point function to first 
order in A. 




Figure 6: The two-point function to first order in A. 

bosonic or fcrmionic, provided the fields satisiy local canonical commutation (or 
anti-commutation) relations. 

The expansion of the two point function has the form 

(</)(xi)0(X2)) = ^ E / '^''y^ ■ ■ ■ '^"yn ^—^{<t>{x^)<t>{x2) n Ant (<A(%)))o 

^ ' ri=0 •' ' j=l 

(37) 

The denominator factor Z[0] has a similar expansion 

= E ^ / '^''y^ ■ ■ ■ ^''yn (R -^-t {Hyj)))o (ss) 

where {A{(j)))o denotes an expectation value of the operator A{0) in the free 
field theory. 

Let us consider first the numerator. Each expectation value involves a sum 
of products of pairwise contractions. If we assign a Feynman diagram to each 
contribution, it is clear that we can classify these terms into two classes: (a) 
linked and (b) unlinked diagrams. A diagram is said to be unlinked if it contains 
a sub-diagram in which a set of internal vertices are linked with each other but 
not to an external vertex. The linked diagrams satisfy the opposite property. 
Since the vacuum diagrams by definition do not contain any external vertices, 
they are unlinked. 

All the expectation values that appear in the numerator can be written as a 
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sum of terms, each of the form of a linked diagram times a vacuum graph, i.e., 

{<i){xi)(j){x2)Ci^t {(t>{yi)) ■ ■ ■ Ant {(t>{yn)))o 

n y \ k n 

= E( jW^l)^(^2)nAntWy,)))^( n AntW%)))0 



fe=0 ^ ' j=l j=fc+l 

(39) 

where the super-index £ denotes a linked factor, i.e., a factor that does not 
contain any vacuum sub-diagram. Thus, the numerator has the form 



n=Ofe=0 ■ ^ ^ j=l j=k+l 

(40) 

which factorizes into 

/ (<^(a;i).^(a;2)n'Cint W%)))o) 

fe=o ■ ^ j=l 

X ( J2 —r / ^""yi • • • ^^''yniU Ant (<^(%)))0 
n=0 j=l 

(41) 

We can clearly recognize that the second factor is exactly equal to the denomi- 
nator Z[0]. 

Hence wc find that we can write the two-point function as a sum of linked 
Feynman diagrams: 

{<l>{xi)cl>{x2)) = ^-^ / d^yi . . . {cP{xi)4.ix2) n Ant icp{yn)))i 

(42) 

This result is known as the linked- cluster theorem. This theorem, which proves 
that the vacuum diagrams cancel exactly out to all orders in perturbation theory, 
is valid for all the A^-point functions (not just for two-point function) and for 
any local theory. It also holds in Minkowski space-time upon the replacement 
(— 1)" <-> i". It holds for all theories with a local canonical structure, relativistic 
or not. 

11.4 Summary of Feynman Rules for theory 
11.4.1 Position Space 

The general rules to construct the diagrams for the A^-point function 
(0|r^(xi) . . . ^(a;jv)|0) in ^'^ theory in Minkowski space and {^{xi) . . . (f){xN)) in 
Euclidean space, in position space are 



k N 
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1. A general graph for the TV-point function has A'^ external points (ore ver- 
tices) and n interaction vertices, where n is the order in perturbation 
theory. Each vertex is a point with a coordinate label and 4 lines (for a 
(f)^ theory) coining out of it. 

2. Draw all topologically distinct graphs by connecting the external points 
and the internal vertices in all possible ways. Discard all graphs which 
contain sub-diagrams not linked to at least one external point. 

3. The following weight is assigned to each graph: 

(a) For every vertex a factor of —i^ in Minkowski space and — in 
Euclidean space . 

(b) For every line connecting a pair of points zi and Z2, a factor of 

(2) 

, Z2) in Minkowski space, or 
{(j){zi)<p{z2))o = Go{zi — Z2) in Euclidean space. 

(c) An overall factor of ^ . 

(d) A multiplicity factor which counts the number of ways in which the 
lives can be joined without altering the topology of the graph. 

(e) Integrate over all internal coordinates. 
For example, the 4-point function 

G'^^\x-i,X2,X3,X4) = {(l>{xi)(l){X2)(f>{X3)(j){x4)) (43) 

has the contribution at order shown in the diagram of figure 7. These two 



Xi X2 X\ X2 




Xi X4 Xi 



Figure 7: Two contributions to the four-point function to order A^. 

diagrams have exactly the same weight (if Gq (1,2) = Gg "^(2,1)), and their 
total contribution to the 4-point function is 

h. (ir)'^ / '^''y' j '^''y^ ^o'H^i,2/i) G^i^\x2,yi) [G'i\yuy2)\^ G'i\xi,y2) G'i\x4,y2) 

(44) 

but are topologically distinct and, thus, count as separate contributions. The 
multipUcity factor S is 5 = (4 x 3)^ x 2. 
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11.4.2 Momentum Space 

The A^-point functions can also be eomputed in momentum space. The rules 
for constructing Feynman diagrams in momentum space arc: 

1. A graph has iV external legs, labelled by a set of external momenta 
ki,...,kN, flowing into the diagram and n internal vertices (for order 
n in perturbation theory). Each vertex has 4 lines (for tp'^ theory) each 
carrying a momentum qi,...q4 (out of the vertex). All lines must be 
connected in pairs. All vacuum terms have to be discarded. 

2. Draw all the topologically different graphs. 

3. Weight of each diagram: 

(a) for each vertex, a factor of S"^ {^f^-^ qj). 

(b) each line carries a momentum and contributes to the weight with 
a factor G^\p) — a^^^'a , in Euclidean space. In Minkowski space it 
becomes —iG\ = — % , . . 

(c) all the numerical factors are the same as in position space. 

(d) we must integrate over all the internal momenta. 

For example, the first order contribution to the two-point function is the tadpole 
diagram shown in Fig. 8. It has the algebraic weight 



Figure 8: First order contribution to the two-point function in momentum space. 



1 




V 



V 




(45) 
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11.5 The two-point function and the self-energy 

To first order in A, and in momentum space, the two point function in Euclidean 

space is 

(46) 

Let us define by the effective or renormalized mass (squared) such that 

(47) 



Again, to first order in A, we can write the equivalent expression 

G^^\p) = + 0{\^) (48) 

P + 2] {2n)^q^ + m^ 

This expression leads us to define is^ to be 

u.^ =m^ + - / - — -^^5 ^ + . . . (49) 

^ 2 J (27r)^ q^ + m? ^ ' 

This equation is equivalent to a sum of a large number of diagrams with higher 
order in A. How do we know that it is consistent? Lot us first note that we have 
summed diagrams of the form shown in figure 9. These diagrams have the very 




Figure 9: The set of all one-particle reducible diagrams of the two-point function 
to leading order in A. 

special feature that it is possible to split the diagram into two sub-diagrams by 
cutting only a single internal line. Momentum conservation requires that the 
momentum of that line be equal to the momentum on the incoming external 
leg. Thus, once again we have two types of diagrams: (a) one-particle reducible 
diagrams (which satisfy the property defined above) and (b) the one-particle 
irreducible graphs which do not. Hence, the total contribution to the two-point 
function is the solution of the equation shown in Fig 10. 
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Figure 10: The Dyson equation for the two-point function. 



Here the thick hnes are the full propagator, the thin hne is the bare prop- 
agator, and the shaded blob represents the irreducible diagrams, i.e., diagrams 
with amputated external legs. We represent the blob by the self-energy operator 
S(|?), shown in Fig. 11 

Thus, the total sum satisfies the Dyson equation 

G(2) (p) = G^') ip) + G^') ip) E(p) G(2) {p) (50) 
The inverse of G^'^\p), r(^)(p), satisfies 

r(2)(p) = G^^\p)-' - S(p) = / + ^2 _ (51) 
To first order in A, S(p) is just the tadpole term 

m = -U 2! 2 + 0{>?) (52) 

which happens to be independent of the external momentum p^. Of course, the 
higher order terms in general will be functions of p^. 

In terms of the renormalized mass /i^, to order one-loop {i.e., 0(A)) we get 

.^ = m'-E(ri=™> + i/^-^ + 0(A^) (53) 

Thus, we conclude that vacuum fluctuations renormalize the mass. However, 
a quick look at Eq.(53) reveals that this is very large renormalization. Indeed, 
fluctuations of all momenta, ranging from long wave-lengths (and low energies) 
with q ~ 0, to short wave-lengths (or high energies) contribute to the mass 
renormalization. In fact, the high-energy fluctuations, with 3> m^, yield the 
largest contributions to Eq.(53), since the mass effectively cuts-off the contri- 
butions in the infrared, IR g ^ 0. Moreover, for all dimensions D > 2 the 
high-energy (or ultraviolet, UV, q 00) contribution is divergent. If we were to 
cutoff the integral at a high- momentum scale A, in general space-time dimension 

— =— + h^+h2K^+ ... 

Figure 11: Feynman diagrams summed by the Dyson equation. 
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D the diagram diverges as A^~^. In particular, the tadpole contribution to the 
mass renormalization is logarithmically divergent for D = 2 (1 + 1) dimensions, 
and quadratically divergent for £) = 4 dimensions. 

Thus, although it is consistent (and quite physical) to regard the leading 
effect of fluctuations as a mass renormalization, they amount to a divergent 
change. The reason for this divergence is that all wavelengths contribute, from 
the IR to the UV. This happens since space-time is continuous: we assumed that 
there is no intrinsic short-distance scale below which local field theory would 
not be valid. 

There is a way to think about this problem. The problem of how to under- 
stand the physics of these singular contributions, indeed of how the continuum 
limit (a theory without cutofi^) of quantum field theory is the central purpose 
of the Renormalization Group (RG). We will study this approach in detail next 
term. Here we will discuss some qualitative features. From the point of view of 
the RG the problem is that the continuum theory {i.e., defining a theory with- 
out a UV cutoff) cannot be done naively. We will see next term that for such 
a procedure to work it is necessary to be able to define the theory in a regime 
in which there is no scale, i.e., in a scale-invariant regime. This requirement 
means that one should look at a regime in which the renormalized mass becomes 
arbitrarily small, /x^ — > 0. As we will see below, this requires to fine tune the 
bare coupling constant and the bare mass to some determined critical values. It 
turns out that, near such a critical point a continuum field theory (without a 
UV cutoff) can be defined. The RG point of view relates the problem of the 
definition of a Quantum Field Theory to that of finding a continuous phase 
transition, a central problem in Statistical Physics. 

However, there are alternative descriptions, such as String Theory, that pos- 
tulate that local field theory is not the correct description at short distances, 

typically near the Planck scale, ^pianck = ^ 10~'^^cm (!), where G is New- 

ton's gravitation constant. From this view point, these singular contributions 
at high energies signal a breakdown of the theory at those scales. 

Before we try to compute S(p), it is worth to mention the Hartree Approxi- 
mation. It consists in summing up all tadpole diagrams (and only the tadpole 
diagrams) to all orders in A. A typical graph is shown in figure 12. 

The sum of all the tadpole diagrams can be done by means of a very simple 
trick. Let us modify the expression for the self-energy to make it self consistent, 
i.e.. 



This formula is equivalent to a Dyson equation in which the internal propagator 
is replaced by the full propagator, as in figure 13. This approximation becomes 
exact for a theory of an A^-component real scalar field 4>a{x) {a = 1, . . . N), with 
0{N) symmetry, and interaction 




(54) 




(55) 
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Figure 12: A typical tree diagram. 



in the large N limit, N ^ oo. Otherwise, the solution of this integral equation 
just yields the leading correction. 



Figure 13: The self-energy in the one-loop (Hartree) approximation. 

Eq.(54) is an integral non-linear equation for So(p). Equations of this type 
are common in many-body physics. For example, the gap equation of the BCS 
theory of superconductivity has a similar form. 

Let us now evaluate the integral in the equation for Eo(p), Eq.(54). Clearly 
So (p) is a correction due to virt,ual fluctuations with momenta ranging from 
zero to infinity. These fluctuations do not obey the mass shell condition p"^ = 
w?. Notice that, at this level of approximation, S(p) is independent of the 
momentum. This is only correct to order one-loop. 

Before computing the integral, let us rewrite Eq.(54) in terms of the effective 
or renormalized mass fj? , 



Clearly, ml is IR divergent for D < 2 and UV divergent for D > 2. Let us now 
express the renormalized mass /U^ in terms of m^, and define Sim? = — m^. 





(56) 



Let us denote by ml the value of the bare mass such that fj,^ = 0: 




(57) 
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We find, 

2 r 2 2 ^ /■ 1 

H = 5m +m^ + - — — p 

2 J (27r)-" + /x^ 

_ . 2^A /-^^ ^ 1 1^ 

- +27 (27r)^ U' + M^ 
A 2_ 2A f d^q 1 

- '"^ 2 7 (27r)^ g2(^2+^2) 

(58) 

In D = 4 dimensions tiie integral contributes with a logarithmic divergence (c./. 
Eq.(92) of Appendix A): 

Thus, although the stronger, quadratic divergence, was absorbed by a renor- 
malization of the mass, a weaker logarithmic singularity remains. It turns out, 
as we will discuss next semester, that this remaining singular contribution can 
absorbed only by a rcnormalization of the coupling constant A. In any case, 
what is clear is that, already at the lowest order in perturbation theory, the 
leading corrections can (and do) yield a larger contribution to the behavior of 
physical quantities than the bare, unperturbed, values, and that this corrections 
are not small. 

We will not give a thorough discussion of these singular contributions right 
now. A complete discussion of this problem involves the development of the 
idea of rcnormalization and of the rcnormalization group, which we will do next 
term. However, it is worth to discuss some of the physical issues behind these 
problems. In a relativistic field theory there is no natural cutoff since a cutoff 
would break Lorentz invariance. However, if a field theory like the present 
one is regarded as an effective theory {i.e., not "fundamental") which is only 
correct at distances larger than some scale ^, we can legitimately cutoff the 
integrals at a momentum A > 1/^. But, in this case, we have to argue that, 
at length scales shorter than ^, there is a consistent theory which is free of this 
divergence. Except from some radical new approach (such as String Theory), 
all theories known to date contain divergencies. Does it mean that they are 
meaningless? For a long time {i.e., from the 30's to the 60's) it was assumed 
that the divergencies signaled that QFT was incomplete. However, in the late 
60's and early 70's a new approach was found that made sense of such theories. 
This approach, known as the Rcnormalization Group, tells us that the apparent 
shortcomings of perturbation theory are due to perturbation theory, not to the 
theory itself. In fact, the Rcnormalization Group gives a framework to define the 
theory at so-called non-trivial fixed points of a certain set of transformations, 
in which the results are physical. From this point of view, the problems are not 
the theories but our clumsy computational tools. 

We will end with a discussion of what these large perturbative corrections 
mean in the context of the theory of phase transitions since, after all, it is 
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also described by a theory with the same form. We saw before that in the 
Landau theory of phase transitions the mass (squared) is related to the difference 
= T — To, the distance to the mean field critical temperature, Tq. We can 
think of ml=Tc — Tq, defining a corrected critical temperature, Tc, 



(27r)^ g2 



(60) 



which shows that fluctuations suppress Tc downwards from Tq due to (precisely) 
the large contributions from short distances. Thus, the UV singular effects can 
be absorbed in a new (and lower) Tc- The subtracted mass, Srri^ = T — Tc, is 
now the new control parameter, as shown in the last line of Eq.(58). 

Returning to the mass renormalization of Eq.(58), we can use the integral 
of Eq.(91) (of Appendix A), to find the result 



/i^ = drn^ - 



r(2 



2 (4^)^/2 D _ ^ 

Since Sm"^ = T — T^ we can write this result in the suggestive form 



Srn^ =T- Tc 



r(2 



2 ) 



2 (47r)^/2 



- 1 



(61) 



(62) 



We will now look at the IR behavior, where the renormalized mass 0. 
Let us recall that the relation between the susceptibility and the effective (or 
renormalized) mass: ~ X~^- For D < 4 the second term of Eq.(62) vanishes 
more slowly than the first term (linear in /x^), whereas for D > 4 the first term 
always wins. Hence, for D < A, the one-loop perturbative correction becomes 
more important than the bare (linear in /i^) term, leading us to expect that 
for Z) < 4 the second term in Eq.(62) should give the important contribution, 
whereas for _D > 4 this term becomes negligible as T — > Tc- Hence, at one-loop 
order, we would predict that 



{T-Tc) "-^ 
X{T)^{ (T-Tc)-' 

(T — Tc)~^ X small logarithmic corrections 



D < 4 

D > 4 (63) 
D = 4 



We see that one effect of these fiuctuations can be to change the dependence of 
a physical quantity, such as the susceptibility, on the control parameter, T — Tc, 
that sets how close the theory is the massless or crziicaZ regime. A key purpose of 
the renormalization group program is the prediction of critical exponents such as 
the one we fund in Eq.(63). We will see that this result is actually the beginning 
of a set of controlled approximations to the exact values. Just as important will 
be the fact that the renormalization group will give a deeper interpretation of 
the meaning of renormalization beyond a process of "hiding infinities" . 



19 



12 The four point function and the effective cou- 
pUng constant 

We will now discuss briefly the perturbative contributions to the four point 
function, 

G^'^\xi,X2,X3,X4) = (0(a;i)(/)(a;2)0(a;3)0(a;4)} (64) 

which is also known as the two-particle Green function. We will discuss its 
connection with the effective (or renormalized) coupling constant. 

To zeroth order in perturbation theory, ©(A"), the four point function fac- 
torizes into a product of all (three) possible two point functions obtained by 
pair-wise contractions of the four fleld operators. 



13 13 1 3 




+0(A) 



2 4 2 4 2 4 



In other words, 

G^^\xi,X2,X3,X4,) = 

Go{xi,X3)Go{x2,X4) + Go(a;i,a;2)Go(a;3,a:4) -|- Go(a;i, a;4)Go(a;2, xa) + 0(A) 

(65) 

As it is apparent, to zeroth order in A, the four point function reduces to just 
products of bare two-point functions and hence nothing new is learned from 
it. We will show next semester that to all orders on perturbation theory the 
four-point function has the following structure: 

G^''^(a;i,a;2,X3,a;4) = 

G(') {xi , X3)G^^^ {X2,X4) + G(') (a;i , X2)G^^'> (xg, X4) + G^^) {xi , X4)G^^^ (2:2, X3) 

+ J d^'yi . ..d^'vi G(2)(a;i,yi)G(2)(a;2,y2)G(2)(x3,2/3)G(2)(x4,2/4) T^^\vuy2,yz 

where the factors of G*-^-* {x, x') represent the exact two-point function, and the 
new four-point function, r*-^-* (yi, ?/2j 2/3; 2/4)7 is known as the four-point vertex 
function. The vertex function is defined as the set of one-particle irreducible 
(IPI) Feynman diagrams, i.e., diagrams that cannot be split in two by cutting 
a single propagator line, with the external lines "amputated" (they are already 
accounted for in the propagator factors). 
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In momentum space, due to momentum conservation at the vertex, F^^^ has 
the form 



rW(pi, . . . ,p4) - (2^)^ (j2p}j r(4)(pi, . . . , 



P4) (67) 



The lowest order contribution to F*-^^ (yi, y2, 2/3, ?/4) appear at order A 

FW(yi,y2,y3,y4)=A + 0(A2) (68) 

depicted by the tree level diagrams: 



r(4) A 
\ / \ / 

\ / \ / 

/ \ / \ 

/ \ / \ 



+0(A2 



which, in momentum space is 

rW(pi,...,p4) = A + 0(A2) (69) 

To one-loop order, O(A^), the four-point vertex function is a sum of (three) 
Feynman diagrams of the form 



P3 
/ 

/ 



\ 

\ 

The total contribution to the vertex function F(4)^ to order one-foop, is 

rW(pi,...,P4) = 

, \^ { [ <Pq 1 .1 ^,,o, 

A \ \ TTT-Tn ?V77 ^^5 ^ + ^"^0 permutations } + 0{X ) 

(70) 

This expression has a fogarithmic UV divergence in D = 4, and more severe 
divergencies for D > A. To address this problem let us proceed by analogy 
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with the mass renormahzation and define the physical or renormalized coupUng 
constant g by the value of r('*)(pi, . . . ,^4) at zero external momenta, pi = . . . = 
Pi = 0. (It is up to us to define it at any momentum scale we wish.) 



g= lim r(4)(pi,...,p4)=r(4)(0,...,0) 



(71) 



This definition is convenient and simple but it is problematic if the renormalized 
mass n"^ vanishes {i.e., in the massless or critical theory). To order one- loop, 
the renormalized coupling constant g is 



f d^q 



1 



0(A3) 



(72) 



Using the same line of argument we used to the dine the self energy, 7'-^-', 
we will now sum all the one loop diagrams, as shown in the figure: 




\ 



\ 



+ 
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This "bubble" sum is a geometric series, and it is equivalent to the replacement 
of Eq.(72) by 

5 = A - 3^ / ,f_l ^ ^ ' + 0{X') (73) 



2 J i2TT)D (^2^^2)2 

or, alternatively, to write the bare coupling constant A in terms of the renor- 
malized coupling g as 



X = g + -^9' 



1 



(2^)^ {q^ + ^i^Y 



0{g') 



(74) 
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where wc have replaced the bare mass with the rcnormahzcd mass This 
amounts to add the insertions of tadpole diagrams in the internal propagators. 
This si a consistent at this order in perturbation theory. 

Written in terms of the rcnormahzcd coupling constant g and of the renor- 
malized mass /x^, the vertex function becomes 



rW(Pi,...,P4)=5-Y / 



(2^) 



1 1 



2 



+ two permutations + 0{g^) 



(75) 



which is UV finite for D < Q. Thus, the renormalization of the coupling constant 
leads to a subtraction of the singular expression for the vertex function. 

After the renormalization of the coupling constant, the singular behavior of 
the integral now appears only in the relation between the bare coupling constant 
A and the renormalized coupling constant g, given in Eq.(74). Clearly there are 
a number of ways to interpret the meaning of this relation. One interpretation 
is to say that at a fixed value of the bare coupling constant A, Eq.(74) relates 
the regulator A (which is a momentum scale) and the renormalized coupling 
constant. In other terms, the effective of renormalized coiipling constant g has 
become a function of a momentum (or energy) scale!. Conversely, we can fix 
the renormalized coupling and ask how do we have to change the bare coupling 
constant A as we send the regulator to infinity. 

It will be useful to work with dimensionless quantities. Since the bare and the 
renormalized coupling constants have units of A*~^, we define the dimensionless 
bare coupling constant u by 

A = A^'^u (76) 



Then 



(77) 



We will vary the momentum scale A and the dimensionless bare coupling con- 
stant u at fixed g. The differential change of the dimensionless coupling constant 
is the renormalization group beta-function, 

Piu) = -^^l (78) 

Hence (for Z) — > 4) 

/3{u) = (4 - D)u - -^u^ + 0(y?) (79) 
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(a) D < 4 



(b) D = 4 



(c) 13 > 4 



Figure 14: The beta-function (top) and the IR RG flows (bottom) for D < A, 
D = 4 and D > 4. The UV RG flows are the reverse. 

The behavior of the Renormalization group beta-function (or Gell-Mann-Low) 
for £) < 4, £) = 4 and D > 4 is shown in Fig.l4. The infrared (IR) RG flows, 
i.e., the flow of the dimcnsionlcss coupling constant u as the momentum scale 
A is decreased, is also shown in Fig. 14 (bottom). The UV flows, i.e., the flow of 
u as the momentum scale is increased, is obtained by reversing the direction of 
the RG flows. 

Clearly for D > A the dimcnsionlcss coupling constant flows to in the in- 
frared (at low energies and long distances). That means that in the IR regime, 
A — > 0, the theory becomes weakly coupled, and perturbation theory becomes re- 
liable in that regime. However, at short distances (or at high energies), A — > oo, 
the opposite happens: the dimcnsionlcss coupling becomes large and perturba- 
tion theory breaks down. Four dimensions is special in that the approach and 
departure from the decoupled limit, m = 0, is very slow, and leads to logarithmic 
corrections to the free field values. 

However, for D < A something new happens: there is a non-trivial fixed point 
at u* where /3(m*) — 0. At the fixed point, the coupling constant does not flow as 
the momentum scale changes. Hence, at a fixed point it is possible to send the 
momentum scale A — > oo and effectively have a theory without a cutoff. Notice 
that even infinitesimally away from the fixed point the UV flows are unstable. 
Also, and for the same reason, at the fixed point it is also possible to go into 
the deep infrared regime and have a theory with a finite coupling constant u* . 
It turns out that this behavior is central for the theory of phase transitions. We 
will come back to the problem of the renormalization group next term where 
we will develop it and discuss its application to different theories. 
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A Integrals 



We introduce a momentum cutoff A and to suppress the contributions at large 
momenta, q ^ A oi integrals of the form 

r.2 



Id 



IJ^\ _ f d'^q ^ ^ 
A2 



where we used a Gaussian cutoff function (or regulator). We will only be inter- 
ested in the regime fj,^ <C A^. Using a Feynman-Schwinger parametrization we 
can write 

roo r -^-a{q'+f,^) 



(2f) 



D 



1 



2 



= e^^^^-L^ / dtt-^e-' 



(47r)^/2 



(81) 



Hence 



where r(z/, z) is the incomplete gamma function, with z = ^ and i/ = 1 — ^, 

T{v,z) = / r-^e"* (83) 



and T{i', 0) = r(z^) is the Gamma function 

/■oo 

r(i/)= / dtf-^-* (84) 

If the regulator A is removed {i.e., if we take the limit A oo), /d(a*^) formally 
becomes: 

^o(M^) = r (l - f ) (85) 

For general r(l — D/2) is a meromorphic function of the complex variable 
D, and has simple poles for i/ = or any negative integer, r(l — D/2) has poles 
for D = 2, 4, 6, . . . 

In = 4 dimensions, v = —1 where T{i') has a pole, the incomplete Gamma 
function at = — 1 is (as ^ ^ 0) 

r(-l,z) = Q -Ini") e-' + j (86) 
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where 7 is the Euler-Mascheroni constant 

7 = - / dte-* Int = 0.5772 . . . (87) 
Jo 

Hence, for <C A^, h is 

/.r^)=^-^tfAv^.^ (88) 



A2 J 167r2 87r2 \^i J IGw^ 

Here we see that the leading singularity is quadratic in the regulator A, with a 
sub-leading logarithmic piece. 

In two dimensions h has instead a logarithmic singularity for /z^ <C A^ 



(89) 



Ay 47r V 'A2 



27r \ u, J An 



A second integral of interest is 



AV J g2(^2_^^2) 

= ^ {id{0) - Id (^!^) ] (90) 



JoilJ' /A'^) is UV finite if D < 4, where it is given by 



In four dimensions, J4, has a logarithmic divergence 



A third useful integral is 

'^^A^>-J (27r)^(.2+ 2)2^ - (9-^) 



(92+^2) 

In the massless limit it becomes 

1 



^^(°)= (4.;^/2(Z.-2)(iP-4) ^""^ 



In four dimensions it becomes 
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